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ABSTRACT 

We have found a statistically significant (99. f %) excess of red (O — E > 2) 
galaxies with photographic magnitudes E < 19.5, O < 21 taken from the 
APM Sky Catalogue around z ~ I radiosources from the Uy catalogue. The 
amplitude, scale and dependence on galaxy colours of the observed overdensity 
are consistent with its being a result of the magnification bias caused by the 
weak gravitational lensing of large scale structures at redshift z ~ 0.2 — 0.4 and 
are hardly explained by other causes, as obscuration by dust. 

Subject headings: Gravitational lensing 

1. Introduction 

Several authors have confirmed the existence of statistically significant associations of 
foreground objects, usually bright galaxies and clusters, with high redshift AGNs (Webster 
et al. 1988, Fugmann 1990, Hammer & Le Fevre 1990, Thomas et al. 1994, Bem'tez et al. 
1995, Seitz & Schneider 1995). These associations have been explained as a consequence of 
the gravitational magnification bias (Narayan & Wallington 1989). This phenomenon arises 
when background sources (like the high redshift AGNs) whose line of sight crosses a matter 
positive fluctuation are amplified by a magnification factor //, what makes them more easy 
to detect in that patch of the sky. Besides, the surface density of background objects behind 
the lens is diluted by a factor The resulting effect is described by the expression 

ua(S) = n' A (S I fi) I fi (Narayan 1989), where ua(S) is the observed cumulative number 
density, // is the magnification, and n' A is the unlensed number density which usually has 
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the form n' A (S') oc (£")", S" being the AGN flux in absence of magnification. This yields 
ua(S) oc fjr a ~ x S a and from this expression it is evident that depending on the slope a, we 
can observe an excess (ua > k-'a) or even a deflect (ua < n' A ) of background sources around 
the lenses. Galaxies and clusters trace the matter fluctuations (up to a bias factor 6), and 
this means that depending on the slope of the number counts of background sources, we 
may observe positive, null or negative statistical associations of these foreground objects 
with background sources like distant AGNs (see Bartelmann 1995, Seitz & Schneider 1995). 

The scales at which the mentioned associations seem to be present vary from a few 
arcseconds (Webster et al. 1988, Thomas et al. 1994) to ~ 30' (Fugmann 1990). This seems 
to indicate that the size of the lenses causing the magnification bias varies from galaxies to 
matter overdensities on scales at least as large as clusters of galaxies or even larger (Pei 
1995, Bartelmann 1995). Therefore, the study of this effect provides information about a 
large range of matter fluctuation scales. 

Radio-selected samples of high-redshift objects present a rather steep slope in their 
number counts-flux distribution. This type of samples is therefore very appropiate for 
trying to detect excesses of foreground objects due to the magnification bias. The optimal 
choice would be using a sample of background objects which is flux limited in two or 
more uncorrelated wavelength bands, what can steepen the number counts-flux relation 
considerably (Borgeest et al. 1991, Bartelmann 1995). In this work we shall study the 
distribution of galaxies taken from the APM Sky Catalogue around a sample of AGNs from 
the Uy catalogue (Stickel et al. 1994). 

2. The Data 

The Uy catalogue of radio-sources is described in Stickel et al. 1994. It has been often 
and succesfully used to prove the existence of correlations between several catalogues of 
foreground objects and distant AGNs (Fugmann 1990, Seitz & Schneider 1995, Bartelmann 
& Schneider 1994). Our initial sample of AGNs is formed by all the QSOs and radiogalaxies 
contained in the catalogue with measured redshifts 0.5 < z < 1.5 and with declination 
6 > 0, what yields 91 AGNs. 

The APM Sky Catalogue is described in Irwin et al. 1994. It provides the user with 
files containing all the objects found by the detection algorithm PISA on the scans of the 
Palomar Sky Survey plates. The files contain, among other data, the classification of the 
detected objects (star/galaxy) and their O and E photographic magnitudes. The O band 
is equivalent to U + B and the E band corresponds to a narrow R (Irwin & McMahon 
1992). We asked for the files containing all the objects within a 90' X 90' box centered on 
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each AGN of the sample. Several fields (8) were not available at the APM Sky Catalogue 
database because of the low galactic latitude of the plates containing them (\b\ < 20°). 
Some of the 90' X 90' fields extend over different plates; besides, the plates have f 0' overlaps, 
so some AGNs can be found on two or three different plates. In order to make the fields 
as homogeneus as possible, each of our fields contains data from one sole plate, namely 
from the one which contains the biggest portion of that field. Using only one plate for each 
field means therefore that those fields which happened to be close to a plate border are 
incomplete. In order to have enough statistics to calculate accurately the average density 
within each field, we discarded several fields (10) because they did not have at least 2/3 of 
the required 5400" X 5400" surface within one plate. 

We finally have 73 fields which contain all the objects classified as galaxies in the 
APM Sky Catalogue in the red or blue band and with O < 19.5 and E < 21. The colours 
are consistent to ±0.2 magnitudes and the magnitude faint end zero-points are consistent 
to ±0.25 magnitudes (Irwin et al. 1994). Our chosen magnitude limits are 0.5 mag 
brighter than the limiting magnitudes on each band quoted by Irwin and collaborators. 
Number counts-magnitude plots show that they can be reasonably considered as a sort of 
completion limits in the sense that fainter than these limits the number counts-magnitude 
function suddenly gets steeper and the spurious detections begin to completely dominate 
the number counts. In order to exclude as many spurious galaxies as possible, we require 
our objects to be detected and classified as galaxies in both the red and blue bands. This 
discards 40 — 50% of the objects at magnitudes below our completion limits. At brighter 
magnitudes (E < 18), most of them are objects classified as galaxies in one band and 
as stars or merged objects in the other. This is not surprising if we remember that the 
classification of the objects in the APM Sky Catalogue is an automated procedure and 
that the pixel resolution is rather low. At fainter magnitudes we are also losing many 
galaxies that due to their colours cannot be detected in both bands. Nevertheless, the 
followed procedure is necessary in order to obtain an enough homogeneus sample. Our final 
galaxy sample is therefore formed by the fields centered on 73 AGNs and contains 35, 000 
objects classified as galaxies in the APM Sky Catalogue in the red and blue band and with 
E < 19.5 and O < 21. 

To test for the existence of systematic gradients on the fields, we divided each field 
into four 45' X 45' boxes and added the number of galaxies found on each of these boxes 
over all the fields, after correcting by the 'useful' surface of each box on each field. The 
differences among the number of objects found in the four boxes agrees well with a y/~N 
poissonian dispersion. Another significant test is to consider the number of objects within 
the central part of all the fields containing half of the total 5400" X 5400" surface (excluding 
a central 1000" X 1000" region which contains most of the excess) and compare it with the 
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number of objects on the outer half, after correcting by the useful surface. The difference is 
again consistent with the poissonian y/~N. These two tests are a strong confirmation of the 
homogeneity of the fields on average. 



3. Main results and discussion 



In order to look for density enhancements on a certain scale /, we shall use the following 
method. We divide the fields into boxes of size I 2 and define the density enhancement on 
each box as 

4 



q ^ ~ j. N f 



E 



/ Sf 



where n{ 3 and s 1 - are the number of galaxies and the useful surface respectively in the box 
ij of the field /, is the total number of galaxies of a field and its total useful surface. 

We estimate the statistical significance of the result in a straightforward and 
distribution-independent way through the empirical probability p } p = Ni ess /N tota i where 
Ni ess is the number of boxes (not including the central one) with a value of q less or equal 
than that of the central box and N to tai is the total number of boxes in the field. This 
estimator is quite conservative since it tends to underestimate the real signification of the 
excess as the maximal signification is always p ma x = {N to tai — l)/N to t a i no matter how far 
the value of the central box is from the average. The empirical distribution of q values is 
asymmetric with a tail toward high values of q. We find that p is always similar or smaller 
than the value of the significance obtained using the empirical rms (which is a factor 
1.2 — 1.4 the one deduced from poissonian statistics, Bem'tez et al. 1995) and assuming 
a gaussian distribution. In any case we think that our galaxy fields are large enough to 
provide meaningful values of p ma x for the considered angular scales. The main strenght of 
this method is thus its simplicity and the fact that it does not make any assumptions on 
the probability distribution of the variable q. 

The results for scales of 500"and 1000"and different colours are presented in table 1. It 
can be seen that if we consider all the galaxies with E < 19.5 and < 21 we do not find a 
very significant excess. However, the value of the density enhancement strongly increases 
for red galaxies with — E > 2 and reaches a value of q = 1.298 ± 0.114 (rms error), 
the highest of all the 500" X 500" boxes in the field, corresponding to an empirical p of 
99.1%. This behaviour agrees well with elementary physical considerations. Background 
objects with redshift z s 0.5 — 1.5 are most efficiently magnified by foreground objects 
with z\ ~ 0.2 — 0.4 (Schneider et al. 1992). Most types of M* galaxies at these redshifts 
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will have colours — E ^ 2 (H = 50, A = 0, = 1. M* and k-corrections are taken from 
Driver et al. 1994, and we have used the relationship O — E 2(B — V) quoted in Irwin 
& McMahon 1992 ). Therefore, selecting galaxies with O — E > 2 means isolating the ones 
most likely to be tracing the potential wells acting as lenses for the AGN of our sample, 
what evidently maximizes the excess, otherwise diluted among the objects at lower redshifts 
which belong to structures not contributing so strongly to the lensing effect. The scale at 
which the density enhancement is more significant, 500", corresponds to 2 — 3 Mpc at 
redshifts of z = 0.2 - 0.4. 

In order to show the dependence of the excess with the angular scale we plot the 
cross-correlation AGN-galaxy 

UM{ , )= -i 

where Nq(9) is the number of galaxies found on the annular bin 9 of the field / and (6) 
is its expectation based on the number of objects on the whole field normalized to the 
useful surface of the bin, for the galaxies with O < 21, E < 19.5 and O — E > 2 (figure 1) 
as a function of the radius in arcsec. The error bars are poissonian v^V- 

The previous results should be confronted with the theoretical expectations of 
galaxy-QSO cross-correlations loqg(0) produced by weak and moderate gravitational 
lensing. Recently Bartelmann (1995) has calculated the expected cross-correlation for 
samples of galaxies with different depths and AGN catalogues simulating the Uy, within 
a CDM cosmogony. The value of the parameter q on a scale of 500 " (q — 1 on a scale / 
coincides with the expression for the cross-correlation amplitude within a radius 1/2, except 
that the former is defined on a square box and the latter on a circular surface) in the 
case of no colour constraint on the APM galaxies (see tablel) agrees marginally with the 
results of Bartelmann 1995. In any case, due to the small number of AGNs considered and 
the rather bright magnitude limit of the APM galaxies the amplitude is not significantly 
different from zero as expected by Bartelmann. However it is noteworthy that, as shown 
above, considering only the reddest galaxies produces a strong enhancement of loag as we 
select the objects tracing the potential wells which more effectively act as lenses for the Uy 
sample. 

We have also tried to measure the effect of imposing a double flux limit by setting a 
cutoff on the optical magnitude of the AGNs. Although it seems to be a tendency towards 
greater excesses for AGNs with brighter optical magnitudes, this constraint reduces the 
sample and lowers the statistical significance, what makes the found results not conclusive. 
Larger samples are needed in order to better establish this effect. 
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The observed overdensity of foreground galaxies around high-redshift AGN seems to 
be a result of the gravitational lensing of the large scale structure and is difficult to explain 
by other causes, as obscuration by dust. In fact, it is easy to see that this last effect should 
produce the opposite of what is observed: the objects forming the excess around high 
redshift AGNs would tend to be slightly bluer than the average. 

This research made use of the APM measurements of the POSS plates, accesible on line 
from IoA at Cambridge. The Uy catalogue was obtained from the Astronomer's Bazaar, 
the CDS service for astronomical catalogues accesible on line from Strasbourg. We thank 
J.L. Sanz and J. Silk for interesting comments and discussions and E. Stengler-Larrea for 
technical support. NB and EMG acknowledge financial support from the Spanish DGICYT, 
project PB92-074I. NB acknowledges a Spanish M.E.C. Ph.D. scholarship. 
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FIGURE CAPTIONS 

Figure 1. The cross-correlation function loag between Uy AGNs with 0.5 < z < 1.5 and 
APM galaxies with E < 19.5, O < 21 and O - E > 2. Error bars are la \Jn{0). 



TABLE 1 

Data for the central box 

O — E 500 arcsec 1000 arcsec 

all 1.109 ± 0.078, 91.4% 1.028 ±0.041, 71.7% 

>1 1.122 ± 0.082, 92.3% 1.039 ± 0.047, 83.9% 

> 2 1.298 ±0.114, 99.1% 1.089 ± 0.059, 91.7% 



Note. — Values of q and p for different scales and colours. The errors in q are the empirical rms over all 
the boxes in the field 
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